Topological confinement in graphene bilayer quantum rings 
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We demonstrate the existence of localized electron and hole states in a ring-shaped potential kink 
in biased bilayer graphene. Within the continuum description, we show that for sharp potential 
steps the Dirac equation describing carrier states close to the K (or K') point of the first Brillouin 
zone can be solved analytically for a circular kink/anti-kink dot. The solutions exhibit interfacial 
states which exhibit Aharonov-Bohm oscillations as functions of the height of the potential step 
and/or the radius of the ring. 

PACS numbers: 73.21.-b, 71.10.Pm, 73.23.-b 



Graphene, a one atom thick crystal sheet of carbon, 
has been shown to display striking electronic and me- 
chanical properties which are expected to lead to the de- 
velopment of new devices (for a review, see e.g. [l|). 
These properties are related to the unusual electronic 
structure of graphene, in which the charge carriers be- 
have as massless fermions with a gapless linear disper- 
sion. For two coupled graphene sheets, known as bilayer 
graphene (BG), the electronic structure is modified due 
to the interlayer interaction, with the otherwise linear 
dispersion becoming approximately parabolic. Another 
important feature of BG is the fact that the electronic 
dispersion can develop a gap, either by doping of one of 
the layers or by the application of an external perpendic- 
ular electric field. Such a gap can be tuned by varying 
the external electric field which allows for the possi- 
bility of tailoring the electronic structure of BG for the 
development of devices, such as quantum dots [1, and 
quantum rings 0| . Recently it has been shown that an- 
other consequence of the existence of a tunable gap in BG 
is the possibility of topological confinement of carriers in 
antisymmetric potential 'kinks', i.e. at the interface be- 
tween two regions of an antisymmetric external electric 
field 0-Q ■ These states have similarities with the surface 
states of topological insulators [13, [HI- Their energies are 
found inside the gap and the wavefunctions are predicted 
to decay away from the interface of the kink potential. 
These topological states are expected to be robust with 
respect to the effect of disorder, with the carrier propa- 
gation along the potential kink displaying electron-hole 
asymmetry. In addition, the geometry of the potential 
interface is determined by the shape of the voltage gates 
used to induce the gap. That in turn can be used to 
further constrain the carrier propagation. 
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In this letter we propose a system in which topological 
confined states are realized at the interface of a kink po- 
tential shaped as a ring. This system can be regarded as 
a good approximation of an ideal quantum ring of zero 
width. We obtain analytical expressions for the electronic 
wavefunction of the BG four-band Hamiltonian. Bilayer 
graphene can be described as two bipartite coupled sheets 
with four triangular sublattices labeled as A {A') and B 
{B') in the upper (lower) layer. For Bernal stacking [l^ 
the coupling between layers is described by the hopping 
energy t — 400 meV between sites A and B', so that the 
Hamiltonian around the K point of the first Brillouin 
zone can be written as 
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where Ui and U2 are external electrostatic potentials ap- 
plied respectively to the upper and lower layers. In po- 
lar coordinates we have tt = —ihvpe^^ (^■§^ + and 

TT^ = —ihvFe~'^^ (^-^ — The eigenstates of the 

Hamiltonian Eq. ([1} are four component pseudo-spinors 

^ ^ [i>A, ipB, ipB', 4' A'] , where ipA and ips {^Pa' 
and V'-B') are the envelope functions for the probability 
amplitudes for sublattices A {A') and B {B'), respec- 
tively, in the upper (lower) layer. The resulting four cou- 
pled differential equations can be decoupled, giving, for 
U2 = -Ui 



V^VVa + 2(e2 + C/2)VVa+ 



(2) 



with Ui = Ui/hvp, U2 — U2/fivp, t — t/hvp and e = 
E/hvp, where E is the energy. 

For the case of a quantum ring with radius R, following 
Ref. Q, we assume a sharp potential kink, so that one 
can define two regions, namely I) < p < i? and II) 
p > R. The upper inset of Fig. [T] illustrates a sketch of 
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the system considered in the present letter. In order to 
simplify the calculations, the potentials Ui (red, dashed- 
dotted) and U2 (blue, solid) are assumed to be piecewise 
constant, defined as V and —V for region I and —V and 
V for region II, respectively, as illustrated in the lower 
inset of Fig. 1. Notice that this quantum ring potential 
is different from the one studied in Ref. [s^, where the 
potential was defined as zero inside a finite width circular 
ring and Ui — —U2 — V otherwise. In order to obtain an 
analytical solution for this problem, one must solve the 
system of differential equations for each region and match 
the eigenfunction at the boundaries. For region I the 
solutions are Bessel functions, which obey the relation 
V'^Fm{ap)e'"''^ = ±a^Fjn{ap)e'"'^, where the function 
Fm{ap) denotes the Bessel functions Jm{ap) or Ym{ap), 
with eigenvalue +a^, or the modified Bessel functions 
Im{ap) or Km{ap), with eigenvalue —a^. The circular 
symmetry of the problem implies that m = 0,±1,±2, .... 
By substituting the solutions one finds a fourth order 
algebraic equation for a, whose solutions are 



a 



Henceforth, we consider only energies within the 
interval where a is complex, such that Fm{ap) ex- 
hibits real and imaginary parts. If one chooses a = 



V^ + E^- y/W^E^ -t^{V^ - E^), four linearly in- 
dependent (LI) solutions are found: 3fJ[Jm(ap)]e""^, 
3[J™(ap)]e*'"^ U[Kra{iap)]e'"''> and 5[K™(iap)]e™^ 
It can be verified that these functions are solutions of 
Eq. although they are not separately eigenfunctions 
of the Laplacian. It is important to point out that, for 
this choice of a, this set of functions is LI only when the 
imaginary part of a is negative. 

The functions ^[J„^{ap)y"^'' and 3[ J„i(ap)]e™'' are 
finite at the origin, but diverge as p — 00 for any value 
of TO. On the other hand, 5ft[/fm(*ap)]e*™^ diverges for 
p — 7- and vanishes when p ^ 00, for any value of to. 
The function 5[ifm(*Q^p)]e""^ has the same behavior as 
^[K„i{iap)]e^™^ for to ^ 0, but for to = it is finite at 
the origin. Based on this we construct the wavefunction 
(p^ for the region I as [l^l 

ij^{p,9) = {An[J„,{ap)] + B^[J„,{ap)]}e""' 



(4) 



Using this expression for ip^, the radial part of the 
other components of for region I are 
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FIG. 1: (Color online) Energy states with angular momen- 
tum index m = 0, ±1, ±5 as a function of the square root 
of the kink/anti-kink potential height (lower inset) for a bi- 
layer graphene ring as schematically shown in the upper inset, 
considering a ring radius i? = 50 nm. Full (empty) symbols 
refer to negative (positive) m. The black solid lines depicts 
the limits ±V between the confined states and the continuum, 
whereas the thin dotted line is i? = 0. 



B3[a7i J,„+i(ap)]}, 



(7) 



E)^ and iJ<{p,e) 



where 71 — a — {V 

</)<(p)e*(™-l)^ ^<,{p,e) - 0<,(p)e'Me, and ^< (p,0) = 
4>A' (/3)e*'-™^^-'^. Similarly, for region II, we choose for the 

refunction V'^ 
go to zero as r 



wavefunction ip^ a linear combination of solutions that 



00, namely, 

^>{p, 9) = {C^[K^{iap)] + i?5[A™(mp)]}e""^ (8) 
and find the other components of as 



{-C'^[aKm^i{iap)] + 



{V + E) 
DlR[aK^-i{iap)\) 



-1 



t{V + E) 
D9[72-ftr,„(ia/9)]}, 



{C^Yt2K^{iap)] 



(9) 



(10) 



V-E' 
1 



tiV-E) 



ip) 



(5) 

{Adi[-fiJ^{ap)] + B^[jiJra{ap)]}, 

(6) 



{A'Si[a-fiJ„^+i{ap)] + 



t(y2 _ E2) 



4>1>{P) 



t{V^ - E^) 



{~C'^[aj2Km+i{iap)] + 
+i(«ap)]}, (11) 



where 72 =a'^-iV + E)^ and V'b(p,6') = </'B(/3)e'("'"^)^ 
<l>%<{p)e'^"'^' and V^,(p,0) = <P^,{p)e^^^+'^' . 
The continuity of the wavefunction at p — R implies that 
^<{R^e) = *>(i?, 6*). That condition leads to a system 
of equations from which we obtain the energy eigenvalues. 
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• * o^°°^°''°°^'^°*^2^x^^*'^^'''^^''''L^^^ Lz — hm/R — ±hky^\ zero energy states appear, hence, 
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proportional to the square root of the potential height 
V for V t. In the bilayer graphene ring problem, an 
analogy can be made between the angular momentum 
Lz = hm/R (see e. g. Eq. (23) of Ref. [6]) and the linear 
momentum hky of the onc-dimcnsional case. Whenever 
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FIG. 2: (Color online) Energy spectrum for m = 0, ±1, ±7 
as a function of the ring radius for bilayer graphene in the 
presence of a circular kink/anti-kink potential of height V — 
20 meV. Full (empty) symbols refer to negative (positive) 
values of m. Three lines are drawn to help visualization: the 
black sohd lines delimits the potential height ±1^, whereas 
the thin dotted line is £ = 0. 



The energies of the topologically confined states of a 
bilayer graphene ring with radius R = 50 nm are shown 
in Fig. [T]as function of the square root of the kink/anti- 
kink potential height. The energy spectrum exhibits 
the symmetry E{m) = —E{—m), which corresponds 
to E{ky) = —E{—ky) as found for a one-dimensional 
kink/anti-kink potential in the si-direction [tJ. Notice 
that the m — state is not the lowest energy state, and 
the less energetic confined electrons in such a system will 
have non-zero angular momentum index, even in the ab- 
sence of a magnetic field. The zero energy states are 
two- fold degenerate (without taking spins into account), 
and are realized only for specific values of the potential 
height V. 

Figure [2] shows the energy states with angular momen- 
tum index m = 0, ±1, ±7 as a function of the ring ra- 
dius R for a kink/anti-kink potential height V = 20 meV. 
The solid lines delimits the -20 meV < i? < 20 meV en- 
ergy spectrum of confined states. Notice that double de- 
generate E = states are observed only for specific values 
of the radius. The energies converge to two values, ±Ec, 
as the ring radius increases, forming two merged bands 
around these energies. The value Ec ~ 13.9 meV is iden- 
tical to the energy of the corresponding one-dimensional 
problem with ky = 0. 

The origin of the zero energy states in Figs. [T] and [2] 
is similar to those found earlier for the one-dimensional 
problem: the energy spectrum for a kink/ anti-kink po- 
tential, as shown in Ref. Q, exhibits E ^ states at 
two values of the linear momentum hky with the same 



if one fixes the potential height V, k'^^ will be a fixed 
value, and for each value of m, there will be a value of R 
that satisfies this condition, leading to a double degen- 
erate zero energy state for this value of the radius. For 
example. Fig. [5] shows the results for V = 20 meV, where 
fc^^' is found to be ~ 0.1 nm^^, consequently, zero energy 
states are observed at i? ~ 10 nm for m = ±1, i? ^ 20 
nm for m = ±2 and so on. Moreover, as the E = 
states satisfy the condition m/R = k^y^ ~ \/W /(2'^/'*) 
for F <C i the equally spaced zero energy states ob- 
served in Fig. [T] occur for m/^/y ^ 2-^/'^RVi, which for 
the parameters of Fig. [T] becomes VV ^ 0.9 meV^^^ for 
TO = ±1, ^/V '-^ 1.8 meV^/^ for to = ±2 and so on. It is 
worth to point out that the condition fcy^-* = \/tV /(2'^/*) 
proposed by Martin et al. Q for the one dimensional 
problem was obtained from the reduced 2x2 Hamilto- 
nian, which is valid only for F <C hence the dependence 
of k^y^ on \/V is no longer guaranteed for large values of 
the kink/ anti-kink potential height. 

The fact that the energy of the lowest energy states 
oscillates as function of R and V , where angular mo- 
mentum transitions take place, suggests the possibil- 
ity of observing persistent currents in each valley K 
or K' induced by the external potential, in analogy to 




modulus, say, hk^'y^ and —hky"' , which were shown to be 
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FIG. 3: (Color online) Angular component of the probabil- 
ity density current as a function of the square root of the 
kink/anti-kink potential height V (top), for a ring radius 
i? = 50 nm, and as a function of R (bottom), when V = 20 
meV. The number on each curve refers to the angular mo- 
mentum index m, and the solid curves are a guide to the 
eyes. 
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Aharonov-Bohm rings, but in the absence of a magnetic 
field. The angular component of the probability den- 
sity current [6| for the lowest energy state is shown as a 
function of the kink/ anti-kink potential V (top) and the 
ring radius R (bottom) in Fig. 3, where current jumps 
are observed when angular momentum transitions occur 
between states with different |m|. We point out that 
Jg{m) = Je{—m) and that the current for K and K' val- 
leys have opposite sign, so that the net current, taking 
into account both valleys, is zero. 

In summary, we demonstrated that confined quantum 
ring states can be realized in a circular kink/anti-kink 
potential in bilayer graphene. We obtained an analytical 
solution for the Dirac equation describing electrons close 
to the Dirac point. Zero energy states, with two-fold 
degeneracy, are realized for special values of the radius 



and potential height. Angular currents for the lowest en- 
ergy state, which present oscillations due to angular mo- 
mentum transitions as V or R increases, are observed. 
Although, for the sake of simplicity, the potential pro- 
file was assumed to be abrupt in the present work, in 
a more realistic description it should have a continuous 
shape. Nevertheless, the present results must give at least 
a good qualitative agreement with a kink/ anti-kink BG 
ring, which would be helpful for the understanding of 
future experiments on such a system. 
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